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Abstract 

A qualitative criterion for a pursuer to intercept a target in a class of differential games is 
obtained in terms of future cones: Topological cones that contain all attainable trajectories of 
target or interceptor originating from an initial position. An interception solution exists after some 
initial time iff the future cone of the target lies within the future cone of the interceptor. The 
solution may be regarded as a kind of Nash equillibrium. This result is applied to two examples: 

1. The game of Two Cars: The future cone condition is shown to be equivalent to conditions for 
interception obtained by Cockaynei^ 

2. Satellite warfare: The future cone for a spacecraft or direct-ascent antisatellite weapon (ASAT) 
maneuvering in a central gravitational field is obtained and is shown to equal that for a spacecraft 
which maneuvers solely by means of a single velocity change at the cone vertex. 

The latter result is illustrated with an analysis of the January 2007 interception of the FengYun- 
IC spacecraft. 
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I. INTRODUCTION 



A variety of pursuit /evasion problems may be treated by the methods of differential game 
theory.-i^ These range from idealized problems of an illustrative nature, such as the Homici- 
dal Chauffeur or Two Cars games,-*^"^ to detailed studies of optimal strategies in air-to-air 
combat." In a previous paper-, a qualitative criterion for interception or capturability was 
devised for a class of differential games of kind suggested by Pontryagin.— The criterion 
is a simple one: Let the future cone be the set of all attainable trajectories available 
to a player subsequent to some initial time. Interception of the target by the interceptor is 
guaranteed to be possible if 

c (V\ 

target ^ interceptor V / 

subject to additional assumptions given below. The condition takes the form of a dominant 
strategy Nash equillibrium, in that, no matter how the target maneuvers, the interceptor 



both. 



based on the Lefschetz 



can always maneuver to intercept it at some point to the future o ' 
Section [TTl provides a new proof of the main result of Ref. [l, 
fixed-point theorem, which relaxes assumptions regarding the convexity of subsets of the 
future cones made in Ref. Q. In Section IIIIl the criterion is applied to the classic problem 
of the game of Two Cars, and to the problem of two spacecraft maneuvering in a central 
gravitational field. 

This latter problem is applicable to the study of a direct-ascent antisatellite weapon 
(ASAT) engaging a spacecraft in low Earth orbit. An account of the January 2007 intercep- 
tion of the FengYun-lC meteorological satellite by a Chinese ASAT appears in Section HVl 
as a worked example. 

Following a discussion of the results in Section |Vl a final Section |VT] presents conclusions. 

II. CONDITIONS FOR GUARANTEED INTERCEPTION 

This paper is concerned with differential games that describe the pursuit of a target with 
position y{t) at time t, by an interceptor whose position at time t is x{t). Both target and 
interceptor can maneuver freely and autonomously, subject to constraints. The evolution of 
X and y is of the form 

| = (2) 
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and 

|^G(.,„). (3) 

where F and G are assumed to be bounded analytic functions, and u = u{t) and v = v{t) 
are piecewise analytic controls. 

A. The Future Cone 

We recall the definition of the future cone of a maneuvering player as given in Ref. Q, 
and review its properties. The net effect of propulsive forces and of external forces such as 
gravity or aerodynamic drag cause the position Y{t) of the interceptor (respectively, target) to 
evolve as it maneuvers. The evolution will be continuous, but not necessarily differentiable. 
In addition, target and interceptor will be subject to physical limitations on their peak 
acceleration or total velocity change Af . Finally, target and interceptor may be expected 
to interact within a compact subset of R^. 

At an initial time to, call the position of the interceptor r(to)- The set of all possible 
histories for r{t) originating at r(to) comprises a topological cone in with vertex r(to)- 
We call the set of points subsequently accessible to the interceptor in the time interval 
(^1,^2)5 with to < ^1 < ^2 the future cone of r(to), written K^iti : t2;r{tQ))M 

Following Ref. [jj, we assume that the subset of either cone lying to the future of its 
vertex is a manifold with compact closure possessing a timelike foliation. We shall denote 
by K^{t;r{to)) a leaf of the foliation corresponding to a time t. In particular, there is no 
assumption that future cones of either target or interceptor, or their leaves, are necessarily 
convex. 



B. Guaranteed Interception 



The properties of optimal solutions to (j2]) and ([3]) leading to interception has been well- 



studied since Ref. 



10; vide. Refs. 



11- 



13l . A sizeable body of literature is devoted to the 



closely related Homicidal Chauffeur and Two Cars games introduced by Isaacs in Ref 



g. Refs. |5|-l9 



e. 



23 



These studies variously treat games of degree or of kind, but all restrict 
the form of the game to facilitate its analysis, by means such as limitation to a linear game, 
or to fixed velocity ratios for the players, or to piecewise constant radii of curvature. 
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We work here with ([2]) and ([3]) in fairly general form, but restrict attention to the more 
limited goal of finding qualitative conditions under which we may be confident that the 
interceptor can force an interception, optimally or no. The winning strategy for the inter- 
ceptor is to choose a trajectory that intercepts the target at some time to the future of to- 
This choice amounts to a mapping $ : — )■ from the set of all trajectories available to 
the interceptor to its desired actual trajectory. A guaranteed intercept will be said to exist 
when an interception solution always exists, no matter how the target maneuvers within its 
future cone. 

C. New proof based on fixed-point theorem for correspondances on nonconvex 
domains 

We use the Eilenberg-Montgomery theoreroi^ii^ to prove a version of Theorem 3.3 from 
Ref. Q giving the conditions for the existence of a guaranteed intercept: 

Theorem 1: Let t^ < t^ < ti < t^. Then, a necessary and sufficient condition for the 
existence of a guaranteed intercept at ti G [to,ti] is 



Proof: Sufficiency: Suppose that, of all the possible trajectories C Ky(ti;y(to)), the 
actual trajectory of the target is y*(t). The mapping from K^ifa '■ t^]x{ta)) into K^ifa : 
tu}',x(ta)) n K:^(tQ : ti;y(to)) = Ky(fo '■ ti;y(to)) is given by the correspondance = 
{y*{t),t G [to,ti]}, which is continuous. Its value is homeomorphic to a one-simplex in R^, 
and is thus acyclic. The cones K^ifa : ti^;x(ta)) and Kyif^ : ti;?/(to)) are differentiable 
manifolds by construction, and thus triangulable.— Each may therefore be regarded as the 
polyhedron of a simply connected simplicial complex and is, in consequence, both an acyclic 
seir^ and an absolute neighborhood retract.— The conditions for the Eilenberg-Montgomery 
theorem are satisfied: A fixed point 



exists for to < ti < ti. Combining this result with the tautological fixed point y*{ti) G 



K{to : t^-y{to)) C K+{ta : C;x(t,)) 



(4) 



x*{t,)e{y*{t),\/te[to,ti]} 



(5) 
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Ky(to : ti;?/(to)) resulting from the target's ability to maneuver freely, we find 

{y*it),\/te[to,h]}\ 

The proof of the necessary condition is identical to that in Ref. Q but is included here to 
make the treatment self-contained. We begin with a condition for the relation between leaves 
of the two cones at a single time. A necessary condition for the existence of a guaranteed 
intercept is that, at the time of intercept ti, 

K^iU;y{h))CKt{U;x{to)). (7) 

for to,ti < ti'. Suppose a guaranteed intercept exists at time tj. Then, every point y in 
Kyiti] yiti)) must coincide with some point x in K^iti] xito)) in order that || x — y ||= for 
at least one pair of values of x and y. Were ([7]) false, there would be some portion of K^iti) 
that lay outside the attainable set of interceptor positions at that time. Thus there would 
be a subset of Ky{ti] y(ti)) for which || x — ?/ ||> 0, Vx G K^itf, x(to)). 

The necessary condition for the entire cone Ky{to : ti;y{ti)) is obtained by transfinite 
induction^^. Let < < tj^ < < ti < ti < and take tj — t as an ordinal. We prove the 
result for Kyitji : ti;y(to)) C Kyit^ : ti;y(to)) and extend to the full set Ky(to : ti;y{to)) 
at the end. 

We begin by showing the necessary condition holds at late times. Suppose that a guar- 
anteed intercept exists at time tj for t^,ti within any neighborhood of ti. As — )■ ti, 

ir+(t,:ti;y(to))^i^+(ti;2/(to)). (8) 

By the condition for a single leaf, it follows that 

y(^o)) C K+{ti; x(t^)) C K+{ta : t^; x(t„)) (9) 

Next, suppose that at least one guaranteed intercept opportunity exists for time ti be- 
tween t^ and ti. By the inductive hypothesis, 

K+{t^ : ti; y{to)) C K+{t^ : t^; x{t^)) (10) 

We wish to examine the prospects at an earlier time t/j. Consider the sets Ky{tj3 : t^] y(to)) 
and K^itp : t^; x{ta)): 

K+{tp : t^-y{to)) = K+{tp : t^-y{to))UK+{t^ : ti;i/(to)) (H) 

5 



and similarly for K^. But if a guaranteed intercept is to be possible Vt G (t/?, t^), at no time 
t in {tp,t^) can it be that 

K^{t-y{t,))<lKt{t-,x{t^)), (12) 

by the single-leaf condition. Letting — )■ to, we have dl]). □ 

As noted in Ref. [ij, the result (E]) may be interpreted as a Nash equillibrium.— The 
interceptor strategy in (|6]) is a strictly dominant one.— That is, while the strategy of the 
target allows it to move anywhere within its future cone, the interceptor always has available 
to it the strategy which places it at some future position of the target. On the other hand, 
no alternative to the target's strategy given in ([H]) will increase its chances of survival. We 
note that, having stipulated in Ky{t; yito)), any such alternative will necessarily constrict 
the volume available for maneuvering by the target, with (we may suppose) the effect of 
worsening its prospects. If we argue in this way, the strategy in ([6]) is preferable to any other 
available to the target, and thus strictly dominant. Both target and interceptor then have 
available to them a strictly dominant, hence optimal, strategy. We may conclude ([6]) is the 
unique Nash equillibrium for this problem. 

III. CALCULATION OF FUTURE CONES FOR SELECTED PROBLEMS 
A. Two Cars 

In this game, introduced by Isaacs^, two cars maneuver in a plane. We have a target Car 
2 with position {x2{t),y2{t)) at time t, pursued by an interceptor Car 1 whose position at 
time t is {xi(t),yi(t)). Players move with constant velocities Vi and V2, respectively, and 
maneuver exclusively by steering. The motion of Car 1 in the plane is given by 

^ = V, s^n{9^{t)) (13) 



and 



^ = vrc0s{er{t)). (14) 



and similarly for Car 2. The control law for ^i(t) is subject to the constraint 

9, < ^. (15) 
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Cockayne^ showed that sufficient and necessary conditions for Car 1 to intercept Car 2 

are 

Vi > V2. (16) 

and 

We show that the condition for guaranteed interception given by applying Theorem 1 

C Kt (18) 

to the game of two cars is equivalent to the paired conditions (fT6|) - (fT7|) found by Cockayne. 
It is convenient to assume that the initial velocity vector of either car may point in any 
direction. This assumption ensures that the future cones are simply connected. 

To see equivalence of the sufficient conditions, we adapt Isaacs'^ construction of trajecto- 
ries for Car 1 that result in interception of Car 2. Amongst the admissible trajectories G 
must be those produced by the Method of the Explicit Policy. The account of the method 
as applied to this problem given by Isaacs^ can hardly be bettered (P is the pursuer, Car 1 
and E, the evader, is Car 2): 

// P has the higher speed and at least as favorable a curvature restriction as E, capture 
can be attained. For P can first go to E's starting point and then follow his track. 

It is clear that, if C , Car 1 can accomplish this for every trajectory of Car 2 
e K^. Let K^^p be the cone comprised of all trajectories constructed by Explicit Policy: 

Kt D Kt,p = Kt. (19) 

The resulting sufficient conditions required for the motion of Car 1 are two^: 
1. Vi > V2: Assume 

Kt{t2 : t/; r2(ti)) C Ktiti : tj; r^{to)). (20) 

with 

to < tl < t2 (21) 

and tf > t2 large, but otherwise arbitrary. In (I2T1) . one may suppose that the difference 
between t2 and ti includes an overestimate of the time for Car 1 to come about, should its 
initital velocity not be parallel to that of the target. 
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Within K^{ti : t/;ri(to)) consider the future cone for Car 1 originating at r2(ti). We 
wish to force an interception for t > t2. By hypothesis, 

Kt{t2 : t- V2{ti)) C Kt{t2 : t; ri(t2)). (22) 

Car I's position at time t may be written formally as 

xi{t) = vi f sin{9i{s))ds (23) 

JQ 



and 



yi{t) = v^ [ cos{9i{s))ds. (24) 

^0 



The distance traveled from the origin at time t is 

n{t)=viy^^^Tf, (25) 



ViMj^stn{e,{s))dsy +(^ j;^cosie,{s))dsy. (26) 



The Schwartz inequality gives 



(^J^ stn{ei{s)) ds^ < stn\ei{s)) ds ds (27) 



We thus find 



ri{t)<Vi^J^ ds^ J^{sin'^{ei{s)) + cos^{ei{s)))ds (28) 

or 

ri{t)<viti. (29) 

By ( p9|) . the leaf i^^^(t2) of the target at a time t > t2 is contained within a circle of 
radius R2, and the leaf K^{t2) within a circle of radius Ri. From f l22|) . 

R2 = V2{t - t2) < Rl = Vi{t - t2) (30) 

or 

■Wl > V2. (31) 

2 2 

2. ^ > The Method of the Explicit Policy requires 

02 < e'l (32) 

From (fT3|) and (fTil), 

^ = t;icos(^^i(t))0i, (33) 
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Vi sin{ei{t))e 



(34) 



and for Car 2, likewise, whence the acceleration 



a = 



iyi cos{9i{t)) OiY + {vi sin{ei(t))eiY 



(35) 
(36) 




(37) 



by (|T5|) . Note that a ■ v = at all times. Equations (|32|) and (|37|) immediately give 

The sufficient conditions fl3T|) and fl38l) are those found by Cockayne. 

We may see the equivalence of the necessary condition for guaranteed interception 



and ffT6|) - (fT7|) by considering the contrapositive: If the conditions be violated, then by 
Cockayne's analysis, there are points the Car 2 can reach that Car 1 cannot (at certain 
constructable times). If that be true, then K2 ^ ■ 

B. Spacecraft maneuvering in a gravitational field 

In this section, we calculate the future cone of a spacecraft maneuvering in a gravitational 
field by impulsive velocity changes. The result yields an account of a direct-ascent ASAT 
engagement with a target satellite. It is also applicable to the endgame of a co-orbital ASAT 
attack.— 

The motion of the spacecraft is given by 



where r is the spacecraft radius vector, v is its velocity, a(t) its acceleration, and /x = 
GNewt.M^. (Because the spacecraft state descriptor includes the velocity, the condition for 
interception is modified in an obvious way.) We make use of an exact solution for the 



(39) 




(40) 
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system (HOj) in the case of ballistic motion in the Earth's gravitational field. The motion of 
a spacecraft in a bound Keplerian orbit is obtained using Lagrange coefficients 



(41) 



with 

sJ^W'"'*! (42) 
\F,{t) G,{t) 

It is not possible to express the coefficients in (H2i) in closed form as functions of time. 
However, it is possible to do so in terms of the corresponding true anomaly /. If /o is the 
true anomaly for which 

'rm))\Jr. I (43) 

v(t(/o)) ; ' 

O^f-fo, (44) 



and 



we have: 



F = 1 - ^ (1 - cos{e)) (45) 
G = sin{e) (46) 



= 1 - ^ (1 - cos{e)) (48) 



o-o(l - cos{e) - sin{e)\ (47) 



In the foregoing, a is the semimajor axis, e the eccentricity, tq is the initial orbital radius at 
true anomaly / = /o, the parameter p = ^^^^^ = a(l — e^), and 



(To = — —. (49) 
V 

A complete description of the spacecraft motion also requires the eccentric and mean 
anomalies. These quantities are related by 



tani^] = jii±4tanff ). (50) 



2 1 -e V 2 
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(Note that ^ and -| always lie in the same quadrant.) The relation between E and the mean 
anomaly M given by Kepler's equation: 



M = E -esin{E). 



(51) 



By use of the mean motion 



n = 




(52) 



(obtained from Kepler's third law) the time t corresponding to true anomaly / is obtained 
from M and the time of pericenter passage r by 



We wish to calculate the future cone K~^(ti : t2;r(to)) of a spacecraft maneuvering in 
the Earth's gravitational field. To that end, we model the motion of an idealized target or 
interceptor with a trajectory consisting of piecewise ballistic orbital segments punctuated 
by shocks in which 



The acceleration during impulsive maneuvers of actual spacecraft will be bounded and con- 
tinuous, even when occurring in times much shorter than any other timescale characteristic 
of the motion. The description of impulsive velocity changes as instantaneous shocks is an 
idealization adopted for computational convenience. However, in the proof of Theorem 3 
below, the acceleration history of the spacecraft is assumed continuous, however violent, in 
order to satisfy requirements of the Gronwall inequality. In order to include the limiting 
case of instantaneous shocks, we may treat the acceleration history as a distribution. The 
transition in fl55|) may then be regarded as the limit of a suitable sequence of good functions, 
and that limit may be taken at the end of other calculations. 

We assume the spacecraft maneuvers by a number (possibly large) of small-impulse 
shocks, and that the total impulse available to it is limited by 




(53) 



The time interval between two points on a ballistic arc is thus 




(54) 




(55) 




(56) 



11 



where ||x|p = x ■ x. The spacecraft trajectory is thus a chain comprised of piecewise 
eUiptical baUistic orbits that is continuous, but not necessarily differentiable at any point 
corresponding to a shock. For all values of the true anomaly between shocks t„ < t < t„+i 
the motion is given by the free-fall Lagrange coefficients, in expressions to be given, as a 
function oi 9 = f — . In any ballistic interval of the spacecraft motion, / is the true 
anomaly of the osculating elliptical orbit that coincides with the spacecraft trajectory (in 
this connection, by "osculating" one should perhaps understand "liplocked" ). 

We require a relation between the true anomaly and time at any point in the course of 
an engagement. A closed- from solution for the time (found by Kepler) may be expressed 
in terms of the eccentric and mean anomalies. Begin with the relations between true and 
eccentric anomalies within the n*^ ballistic arc for three points on the orbit: 



to, I = ,/(i±£)(„„ (£^ I . (57) 



and 



-(^1-J^-(^1. (58) 



In the foregoing, 9 = f — (respectively, E — Eq) is the difference in true (respectively, 
eccentric) anomalies between the position r and the semimajor axis of the n*^ ballistic arc, 
while 9^ = — J'q {En — Eq) is the difference in true (eccentric) anomalies between the 
position r„ and semimajor axis. The time interval between r and is given by fl54p as 




t-tn = J- {{E- E^) - {E^ - E^q) - e {stn{E - E^) - sin{En - E^^))) . (59) 



At a time tn+i, let a velocity change Av„+i be imparted when the spacecraft is at position 
Vn- The initial position and velocity are thus r„ and v + Av„+i. The semimajor axis of the 
new ballistic arc is obtained from the vis-viva equation 

II V + Av„+if = (-^ - . (60) 

V||r„|| a„+ij 

The parameter (and hence the eccentricity) is given by 

n ||r„ X (v + Av„+i)f 
p = a{l-e) = , (61) 

1^ 



and 



rn • (v + Av„+i) 
(To = — • [o2) 
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These quantities determine the Lagrange coefficient matrix for the n + 1^* balhstic arc. 

Call K^iiti : t2]r{to)) the set of all trajectories originating from r(to) and subject to n 
shocks {Avj, i = 1 . . .n,} respecting fl56p . with the understanding that for n = 1 the shock 
occurs at time to- The future cone of a spacecraft maneuvering by ballistic motion between 
shocks is then 

oo 

K+{h ■■ r(to)) = U K+{h ■■ t2; r(to)). (63) 

n=l 

We now prove a result relating the future cone K^(ti : r(to)) of a spacecraft that ma- 
neuvers in a gravitational field by free fall punctuated with a countable sequence of shocks 
{Avj}, subject to the overall upper limit on their sum ( 156|) . to the set K^(ti : t2;r(to)) 
comprised of the union of single ballistic arcs with ||Av|| < Avtot- If we stipulate that 
the spacecraft be capable of maneuvers for which ||Av|| = Avtot at the cone vertex, a 
particularly simple relation holds: 

Theorem 2: K+{t^ : r(to)) = K+{t^ : r(to)) 
Proof: C K^: By their respective definitions. 

C Kf: By induction. We first show that, corresponding to any piecewise ballistic 
trajectory with n shocks Avj, z = ■ ■ ■ n — 1, an equivalent single ballistic trajectory exists 
connecting the cone vertex with any point on the piecewise trajectory subsequent to the 
final shock. We then show the velocity change Avq at the vertex of the single ballistic arc 
obeys 

n-l 

IIAv^ll < ^ ||Av,|| < Avtot. (64) 

In the case tt, = 1 of a single shock at time t^, it is immediate that 

IIAvoll = ||Av°|| < At;toi. (65) 

Assume that after n shocks the spacecraft location r„(t„) is connected to the vertex by a 
single ballistic arc with '^^Zq l|Avj|| < Avtot- We show that after a further shock obeying 
^"^Q ||Avi|| < Avtot, any point on the resulting ballistic orbit of the spacecraft r„+i(t) with 
t > tn+i is likewise connected to the cone vertex by a single ballistic arc G K^(ti : t2] r(to))- 
Take any point r„+i(t) in (tn+i,t2)- We seek a single ballistic arc connecting the cone 
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vertex and r„+i(t), v„_(.i(t) of the form 



I „ , , I (66) 



with 

Vo"+i=vo + Avo"+i (67) 

for Lagrangian coefficients corresponding to some balhstic orbit C K^{ti : t2;i"(^o))- 

To this relation we append the condition that determines the true anomaly corresponding 
to the final time of the single ballistic arc, 

n+l n 

t{0) -to = $^(ti - tj-i) = [tn+i - K) + XI (^i - h-i)- (68) 
i=i i=i 

Equations ( 166|) and (168|) comprise (in light of the dependence of the Lagrange coefficients 

upon Vq"^^) a nonlinear vector equation in four unknowns of the form 

C(^,v^+\t,r„+i,v„+i) = 0. (69) 

The relation between (6',Vq^^) and (t, r„_|_i, v„+i) is, with stipulated fq, diffeomorphic. 
As a result, a version of the implicit function theorem proved by Kumagai^^ supplies a (real) 
solution of ( 169|) for (6',Vq^^). Thus, there exists a single ballistic arc corresponding to the 
n + l shock piecewise ballistic trajectory. 
It remains to estimate the magnitude of 

Av^+i = Vo"+i - vo. (70) 

The inductive hypothesis is that a single ballistic arc exists corresponding to a trajectory 
which has experienced n shocks, for which 

n-l 

X IIAvill < Avtof (71) 

i=0 

and 

n-l 

||Avo1|<Xl|Av,||. (72) 

Assume that f l72p holds for n > 1. Now consider the full trajectory with n + l shocks, 
subject to the overall limit 

n 

Xl|Av,|| < Ai;toi. (73) 
14 



We may replace that portion of the trajectory corresponding to the first n shocks by the 
equivalent single ballistic arc. The entire trajectory may thus be replaced by one comprised 
of two ballistic arcs. Applying and ( 172]) . we have 

W^^o^^W < W^^oW + ||Av„|| (74) 

n-l 

< J]||Avi|| + ||Av„|| (75) 

1=0 

We conclude 

n 

||Avo"+i<5^||Av,||<At;,„,. (76) 

i=0 

Thus, a single ballistic arc C K^{ti : t2; r(io)) connects each r„ G K:^{ti : r(to)) with 
the cone vertex. Call the set of all such single ballistic arcs K^*. Therefore, 

oo 

K+{h : h; r(to)) = |J K+{h : h; r(to)) (77) 



n=l 



C K+* C Kt, (78) 



proving the theorem. □ 

is a subset of B{[tQ,ti]), the set of bounded functions on [toj^i], which is a complete 
metric space in the sup norm. is also closed: The mapping from 



(79) 




to 

' r„(t) 



(80) 



v„(t) 

in = is homeomorphic (diffeomorphic, in fact) and thus maps closed sets onto closed 
sets. The set {vq} that generates is delimited by Avg < ||Ai;tof||; it, and thus K^, is a 
closed set. The latter is therefore a complete metric space. 

Theorem 2 gives us the future cone for a countable number of shocks in terms 
of the trajectories resulting from a single shock experienced at rito) = tq. The generic 
multishock trajectory amounts to a Devil's staircase in Av proceeding from Av = to 
Av = Av final < Avtot- The theorem to be proved next extends this result to the case of a 
spacecraft maneuvering by continuous nongravitational acceleration as, for instance, in the 
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case of the powered flight of a rocket: 



Theorem 3: Let r{t),t G [toji^i] the trajectory of a continuously accelerating spacecraft, 
originating at yq with initial velocity vq, and subject to Av < Avtot- Then r{t) G K^{ti : 
t2;r(to)) 

Proof: Let r{t) be the trajectory, originating at r(to) with velocity vq, that results from 
an imposed acceleration 

F 

a=— . (81) 

m 

We assume that a is a integrable function of t. The motion of the spacecraft is given by 
TQl): 



d I r \ I V 

di 



+ a(t) 



(82) 



The effect of continuous acceleration is approximated by a series of shocks Avj from 
which we construct a minimizing sequence of shock histories.— Let the trajectory resulting 
from the nth choice of shock history be given by 



d I rn \ I 

dt 



-/i^f + a„(t) 



(83) 



with SLn{t) a collection of impulsive velocity changes punctuating ballistic motion subject to 
gravitation alone. We may approximate a„()f:) for t G [to^i] by a sequence of step functions 

i=l 

as fc„ — > OO, where Xa{x) is the characteristic function for an a subset x G A and the support 
of ajj is suppi^al^) = Ti. The acceleration a„(t) integrates to Av(t) of the form 

/ dsan = J^/i(ri)a^ (85) 

where /i(x) is the (upper) measure of the set x C [toji'^i]- Writing 

r = r„ + (5r 

V = v„ + 5v (86) 
a = a„ + (5a 



16 



we have 



(87) 



Expanding to second order in — 



St 



dt \ / I -ii{5v - y-^v) + 5a(t) / I F(t) 



with 



^/ N A* / r • 5r ^ 3 5r • 5r 15 (r ■ (5r)^ \ ^ f ,5r^..\ . , 

Fm = -4 3^^5r + r + O — M (89) 

V ^ r3 V r2 2 r • r 2 (r • r)2 y / 

We note that control laws for the shock history Avj exist that drive 6r to small values 

6r 

— < 1. (90) 

r 

It may thus be assumed without loss of generality that if 

||5r|U<a||r|U (91) 

in the sup norm on [to;^i] then 3a > such that the second-order terms in flHHl) dominate 
the error resulting from linearization. Using ( |9T|) and the Schwartz inequality, we have on 

Sr ■ Sr < 3\\6r ■ 5r||^ 

II IIOO ^g2) 

and 



r ■ (5r < Sallr ■ r 



|2 

I oo" 



We also have 



6r-6r <^/3\\6r\\^Y.^^ri ^^^^ 
< y/3a\\r\\ooY.i^^i- 



Then, using the estimates (]92il and (l93i) . a reliable overestimate of the second-order error 
that is linear in \6r\ is 
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Upon substitution from and replacement of 



by 



f l88|l takes the form 



r • 6r 

3 — 



gi"- \Sr\ 
3 — ^r, 



(95) 



(96) 



5v " Uv 




(97) 



which, upon formal integration from to to t becomes 

A ' 

where elements of the matrix A(t) > on [to,ti] and we are at liberty to choose a„ so that 




< I dsA{s) ' 

'to 



(98) 



6w{t) 



ds 5a{t). 



(99) 



to 



is positive on [toj^i] for sufficiently large n. 

Using r„ = r + 5p in place of r = r„ + 5r in (!86|) . we find the expected changes to (!88|) : The 
terms linear in 6p change sign compared to (l88i) . while the form of F{t) remains unaltered. 
Following the same development leading to (p8l) . with the change that 



now replaces 



and setting Sp = —Sr, we obtain 



3^r 



3^r, 



(100) 



(101) 



(102) 



with the identical form for the positive matrix function A(s). Both 6r and —6r obey the 
same vector inequality. It must be the case that \6r\ also obeys the inequality (!98|) . The 
multivariate version of the Gronwall inequality — then gives 





"* / 

< / dsV{s,t) 

'to 



6w(s 




(103) 
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where V(s,t) satisfies 



In particular, fllOSp gives us 



V(s, t) = l + A(r)V(r, s) dr. 



(104) 



\Moo < 







f ^ 




< 


j' dsY{s) 1 




1 






[6w I 





(105) 



An elementary but ugly calculation using the Schwartz inequality and the mean value the- 
orem for integrals converts (llOSp into the Lipschitz condition 



M^<c\\Sw\ 



(106) 



Examine now ||(5w||oo- We have, from fl85p 



Sw{t) = lim y /i(rj)a^ — / dsa.{s). 



(107) 



But, a being integrable, we may choose a shock history approximated by a sequence of step 
functions such that^ 

" 108) 



kn nf 

lim V/i(ri)a^ ds . 



Therefore, m exists such that for any 6 > such that 



there is an e > such that 



|5w„||oo < 5, > m. 



||(^r||oo < Vn > m. 



(109) 



(110) 



The set X = {6r,Wn > m} is a bounded subset of the set C([to,ii]) of continuous 
functions defined on [toj^i]- By f lllOp . the error \6r\ is uniformly bounded on [to,^i]- Thus, 
if we choose a neighborhood Nt of any t G [to, ^i], for all t E Nt 



\5y{t) - 5Y{t)\ < \5y{t)\ + |5r(t)| < e 



111) 



uniformly on [to,^!] for all 5y E X. X is therefore an equicontinuous set. The Arzela-Ascoli 
Lemma then implies that a uniformly convergent subsequence (5r„j, exists on [^0,^1] whose 
sup norm tends to zero.— Therefore, as Uk — )■ c>o, r„^.(t) — )• r(t) uniformly on [to,ti]. Recall 
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that is a complete metric space in the sup norm, and thus contains the hmit of all its 
convergent sequences. We may now claim 

r(t) e (112) 

proving the theorem. □ 

Clearly, the future cone of a maneuvering spacecraft is the union of piecewise ballistic 
trajectories and continuously accelerating ones, which gives us 



Theorem 4: K+{ti : r(to)) = K+{ti : ^2; r(to)) 
Proof: C K^: By their respective definitions. 

C K^: By Theorems 2 and 3. If K^^^ is the set K 
{continuously accelerating r(t) G K^}, U K^^^^^ = C K^. □ 
With Theorem 1, Theorem 4 gives immediately a result long known to duck-hunters. 



IV. APPLICATION TO FENGYUN-IC INTERCEPTION 



On or about 22:26 UTC 11 January 2007, a direct-ascent antisatellite weapon launched 
from the territory of the People's Republic of China intercepted, and destroyed, the Chinese 
FengYun-lC weather satellite.— FengYun-lC was launched from the Taiyuan Satellite 
Launch Center on 10 May 1999 into a Sun-synchronous 860 km, 98.8° orbit. It was designed 
for a two-year mission duration. A replacement spacecraft, FengYun-lD, was launched in 
May 2002, but FengYun-lC remained operational into 2005 and evidently responded to 
ground commands as late as January 2007.— 

While the instrument suite and data products of the FengYun meteorological instruments 
are documented in numerous sources, not many descriptions of the spacecraft have appeared. 
FengYun-lC is variously described as having a mass between 880 and 960 kg. The FengYun- 
1 series spacecraft were three-axis stabilized by a combination of reaction wheels and cold 
nitrogen jets. A 1999 report from the Foreign Broadcast Information Service cites an account 
of FengYun-lA appearing in the Journal of Chinese Society of Astronautics that claims the 
attitude control system carried 12 kg of nitrogen at launch.— In an account of the debris 
cloud resulting from the destruction of FengYun-lC, Ref. |3J remarks that at launch, the 



spacecraft mass was 958 kg, but at its end-of-life, the mass was 880 kg. One may interpret 
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this as an upper limit for the reaction mass available for attitude control of 78 kg, again 
presumably nitrogen gas.— 

The application of the results from the preceding section to analysis of the FengYun-lC 
interception of January 2007 is straightforward. One estimates the future cone of the target 
subsequent to some fiducial time, granted credible assumptions about its residual Av, and 
compares this to a model of the future cone of the interceptor for a stipulated engagement 
scenario. A practical method for computational purposes is supplied by Monte Carlo. By 
virtue of Theorem 4, it suffices to pick a number of elements of for each player and (so 
to speak) observe the fall of shot. 

The initial position and velocity of FengYun-lC at the start of the engagement is calcu- 
lated from a two-line element set obtained from Spacetrack^^, using the Aerospace Corpora- 
tion's Satellite Orbital Analysis Program (SOAP)^"^ to propogate the FengYun-lC state to 
the target cone vertex. We use a model of the ASAT developed by G. Forden of MIT.— 
Forden's model assumes that the Chinese direct-ascent ASAT is based on the two-stage 
solid-fueled DF-21 IRBM with a small interceptor of mass 600 kg as a third stage. The 
DF-21 stages are modeled as 1.4 m in diameter, with stage 1 being 5.1 m in length, stage 2 
1.8 m in length. A center-perforated grain consisting of a double base propellant is assumed 
for both stages, with an Igp of 225 seconds for the first stage and an Igp of 230 seconds for 
the second. Each stage burns for 36 seconds. The model neglects Av of the third stage. 

The vertex for the ASAT cone was arbitrarily chosen as the point in its fiyout at which 
the ASAT altitude exceeds 104 km, using the same launch azimuth and initial pitch assumed 
in Forden's analysis. At this point in the fiight profile, the missile is late in its second stage 
burn. The total Av used to calculate K^^^j, is obtained from f lll3p using the partially 
expended second plus third stage mass for rrii and the second stage tare plus third stage 
mass for ruf. This procedure presumably overestimates slightly the Av likely available to 
the actual ASAT very near the cone vertex, but should be accurate enough for illustrative 
purposes. 

The change in rocket velocity Av resulting from the exhaust of reaction mass is given by 
the rocket equation 



where Igp is the specific impulse in seconds, go is the gravitational acceleration at sea level, 
and rrii and rrif are the initial and final masses of the rocket, respectively. The ideal specific 




(113) 
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impulse for a cold nitrogen gas jet is given as 76 seconds in Ref. |41|, p. 229. The range 
12 kg-78 kg for the estimated available reaction mass equates to a total Av in the range 
11m/ s to 63m/s, assuming a tare mass of 880 kg for FengYun-lC. The smaller value of 
11 m/ s is used to calculate the FengYun-lC future cone K~^yic ground that this value 

corresponds to an overestimate of the likely Av available to the spacecraft. Even were the 
total Af = 63 m/s initially, by January 2007 the reaction mass would almost certainly have 
been quite depleted. The vertex of the FengYun-lC cone is chosen at the time of ASAT 
launch. 

The superposition of the future cone K\gjrp of the ASAT and that part of the cone 
KpYic FengYun-lC lying within K\gjij. is presented in Figures [T] and |2] as random-dot 
stereograms with approximately antipodal viewing geometries. The nominal encounter time 
22:26:00 occurs at 450 seconds TALO. In the interval during which the cones intersect, 

Kty.c C K\sAT (114) 

and Theorem 1 guarantees the ASAT can intercept FengYun-lC. 

That the ASAT could intercept FengYun-lC is hardly a novel conclusion-this much was 
demonstrated beyond dispute on 11 January 2007. Admittedly, FengYun-lC appears to have 
served as a passive target.— Consider, however, a hypothetical scenario in which FengYun-lC 
maneuvers during the engagement. The disparity in expansion between the cones Kpy-ic 
and K\gjirp underscores the inability of FengYun-lC to evade the ASAT. If anything, this 
exercise understates the vulnerability of FengYun-lC on that date: The cone KpyiQ is com- 
puted assuming that the total Av remaining in the attitude control system would actually be 
available to the spacecraft for maneuver. It is unlikely this assumption holds in practice. At 
a minimum, the actual Av is limited by the total thrust the cold nitrogen jets are capable of 
producing if operated in a configuration that produces net impulse, even on the assumption 
that operators on the ground issue evasive commands. 



V. DISCUSSION 

It was noted in Ref. Q that ([6]) in Theorem 1 amounts to establishing a (dominant 
strategy) Nash equillibrium: The target may move to any position available in its future, 
but so long as the interceptor's future contains that of the target, the interceptor can always 
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maneuver to some future position r(tj) of the target. Cockayne^ uses the telhng phrase 
"against all opposition" to describe the corresponding situation in the game of Two Cars 

n 

{vide. Isaacs Ref. 3, p. 202, as well). 

Theorem 1 was proved for simply connected future cones, but can be applied to certain 
multiply connected ones. If multiply connected future cones of target and interceptor are 
manifolds with a timelike foliation and are the finite union of polyhedra of simply connected 
simplicial complices, one may invoke Theorem 1 for individual polyhedra. Thus, if 

Ki, {to : tN, 4) = ul'.'Ki.JU., : U; rj) (115) 

where each K^^^ • is a simply connected polyhedron and a subset K^j.g (ta : r?) of K^^g 
is a simply connected polyhedron (or other acyclic absolute neighborhood retract) such that 

KLg.ita : r°) C Kl, Jt,_^ : U,; rj) (116) 

for some j G 1,A^ then the existence of a guaranteed interception opportunity follows in 
any polyhedron of the interceptor cone that contains K^arg.iia '■ C;i"?)- This approach is 
presumably necessary for a treatment of the Game of Two Cars that relaxes the simplifying 
conditions regarding initial velocity and cone vertex times made in Section IIII A[ 

Theorems 2-4 need not give a good account of the future cone of a spacecraft for which 
the assumption that it may expend the entirety of its available Af in a single maneuver at 
the cone vertex is a poor approximation. In that event C K^{^Vtot) remains true, but 
it may happen that ^ Ki{Avavaii)- A proper treatment of this case with the results of 
Section IIII Bl requires considering a sequence of future cones, for each of which the /S.Vavaii. 
available in a given time interval of its motion is used in place of /S.Vtot- 

VI. CONCLUSION 

The sample problems worked in Section IIIII show that the future cone construction 
based on Theorem 1 offers a simple method of determining when a class of differential 
pursuit /evasion games of kind, in which both players maneuver freely, is guaranteed to have 
interception of the target possible outcome. 

The method of analysis presented in Section IIII Bl is applicable to direct-ascent and co- 
orbital ASAT engagements. The specific example of the FengYun-lC interception by a 
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direct-ascent ASAT demonstrates the applicability of the method to the study of a pur- 
suit/evasion game of considerable practical interest. The formalism developed in this paper 
and its predecessor may also be taken to apply, if only as a reducto ad absurdum, to directed 
energy weapon attacks with a free line-of-sight to the target satellite. 
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FIG. 1: Future cones Kpyici'^'^^ • ^'''^i ^(0)) ^^'^ ^asat(^^ • '^^^'^ r(68)) for target and interceptor, 
respectively, in an Earth-centered inertial frame with times in seconds after interceptor launch. 
The stereogram also shows the latitude and longitude lines of the assumed interceptor launch from 
28.13°Af, 102.02°^ (launch azimuth 345.73°) at 22:18:30 UTC 11 January 2007 and the interceptor 
trajectory from launch point to K^sAT vertex. The typical interceptor cone trajectory points into 
the page, so that the launch point appears nearer than the interception region. Distal termini for 
Ktc.,rr are truncated at altitudes below 90 km. 




FIG. 2: Future cones KpY_i(j and K^g^j, for target and interceptor, respectively, as in Figured! 
In this view, the typical interceptor cone trajectory points somewhat out of the page. 
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